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The  problem  under  discussion  is  that  of  calculating  magnetic  field  configurations  in  which  the 
Lorentz  force  j  x  B  is  everywhere  zero,  subject  to  specified  boundary  conditions.  We  choose  to 


represent  the  magnetic  field  in  terms  of  Clebsch  variables  in  the  form  B  =  Vo  x  Vfi.  These 

i***  -  * 

variables  are  cc^ns^ant  on  any  field  line  so  that  each  field  line  is  labeled  by  the  corresponding  values 


of  a  and  0.  When  the  field  is  described  in  this  way,  the  most  appropriate  choice  of  boundary 
conditions  is  to  specify  the  values  of  a  and  0  on  the  bounding  surface.  We  show  that  such  field 


configurations  may  be  calculated  by  a ’’ magneto-frictional”  method.  We  imagine  that  the  field  lines 
move  through  a  stationary  medium,  and  that  each  element  of  magnetic  field  is  subject  to  a  frictional 
force  parallel  to  and  opposing  the  velocity  of  the  field  line.  This  concept  leads  to  an  iteration 
procedure  for  modifying  the  variables  or  and  0,  that  tends  asymptotically  towards  the  force-free 
state.  We  apply  the  method  first  to  a  simple  problem  in  two  rectangular  dimensions,  and  then  to  a 

J 

problem  of  cylindrical  symmetry  that  was  previously  discussed  by  Barnes  and  Sturrock  (1972).  In 
one  important  respect,  our  new  results  differ  from  the  earlier  results  of  Barnes  and  Sturrock,  and 
we  conclude  that  the  earlier  article  was  in  error. 


I.  INTRODUCTION 


Magnetic  fields  play  a  key  role  in  all  solar  activity  and  probably  also  in  similar  phenomena  on 
other  stars  and  in  other  astrophysical  objects.  Knowledge  of  the  quiescent  structure  of  the  magnetic 
field  is  obviously  an  indispensable  prerequisite  for  any  study  of  the  dynamics  of  these  fields.  In 
general,  the  equilibrium  configuration  involves  a  balance  of  several  forces  including  the  Lorentz 
force,  the  plasma-pressure  gradient,  and  the  gravitational  force  acting  on  the  plasma.  However,  in 
many  cases  we  are  dealing  with  a  "low-beta”  situation  in  which  the  gas  pressure  is  small  compared 
with  the  magnetic  pressure.  Since  the  distribution  of  gas  along  a  magnetic  field  line  is  determined 
by  hydrostatic  considerations,  the  gravitational  potential  energy  density  will  be  comparable  with 
the  plasma  pressure.  Under  these  conditions  the  magnetic-field  configuration  will  approximate  a 
force-free  state  for  which  the  Lorentz  force  f  is  zero,  i.e., 

f  =  jxS  =  0.  (1.1) 


It  is  clear  that  the  general  problem  represented  by  equation  (I. I)  is  nonlinear,  since  the  sum  of 
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two  separate  solutions  of  equation  (1.1)  does  not,  in  general,  produce  a  third  solution.  However,  a 
special  class  of  force-free  fields  may  be  generated  from  solutions  of  the  linear  equation 

Vx8  =  \8,  (1.2) 

where  A  is  taken  to  be  a  constant  (Priest  1982).  In  this  case,  it  is  clear  that  the  sum  of  two  solutions  of 
equation  (1.2)  generates  a  third  solution.  As  a  result,  one  can  develop  a  set  of  solutions  of  equation 
(1.1)  in  terms  of  convenient  functions,  and  then  manipulate  a  combination  of  these  solutions  to 
match  specified  boundary  conditions. 

Suppose,  however,  that  we  are  trying  to  find  the  force-free  field  configuration  produced  by  a 
finite  source  such  as  a  pair  of  sunspots  of  opposite  polarity.  If  we  investigate  the  field  at  a  large 
distance  from  the  source,  we  expect  that  the  gradients  will  vary  as  Br~l.  However,  such  a  variation 
of  j  with  respect  to  r  is  clearly  incompatible  with  the  assumption  that  A  is  constant.  Hence  the 
assumption  that  A  =  constant,  which  reduces  the  nonlinear  problem  of  (1.1)  to  the  linear  problem 
described  by  (1.2),  is  inappropriate  for  the  investigation  of  a  wide  range  of  interesting  problems  of 
solar  physics  and  astrophysics. 

In  the  special  case  that  the  magnetic-field  configuration  is  independent  of  one  of  the  spatial 
coordinates,  it  is  possible  to  generate  solutions  of  the  nonlinear  equation  (1.1)  by  a  special  technique 
described  by  Priest  (1982).  For  instance,  for  field  configurations  that  are  independent  of  z  in 
rectangular  cartesian  coordinates,  it  is  possible  to  generate  a  set  of  nonlinear  force-free  fields  by 
specifying  that  the  field  is  of  the  form 

=  =  Bz=F(Az),  (1.3) 

where  F(AX)  is  a  specified  function  of  A  ,  and  the  vector  potential  satisfies  the  equation 

^r  +  ^r  +  fWm)-0-  (i-<) 

However,  this  is  simply  a  technique  for  generating  a  special  set  of  solutions  of  the  nonlinear  equation, 
and  is  therefore  inappropriate  for  solving  a  general  problem,  such  as  computing  the  force-free  field 
that  satisfies  certain  boundary  conditions. 

Force-free  fields  arise  in  situations  in  which  the  electrical  conductivity  is  sufficiently  high  that 
the  relative  diffusion  of  plasma  and  magnetic  field  may  be  neglected.  In  this  situation,  Alfv^n’s 
concept  of  "frozen-in  magnetic  field”  is  appropriate  (Cowling  1957),  and  it  is  convenient  to  think  of 
magnetic  field  lines  as  having  physical  reality.  In  view  of  the  equation 

V-B  =  0,  (1.5) 


it  is  possible  to  describe  the  magnetic  field  in  terms  of  the  Clebsch  variables 


from  which  it  is  obvious  that 

&  •  Va  =  0,  S-V0  =  0,  (1.7) 

so  that  a  and  0  may  be  used  to  ’’label”  the  field  lines  of  the  configuration  (Parker  1979). 

It  is  worth  noting  that  one  could  employ,  instead  of  equation  (1.6),  the  equation 

g=f(a,0)VaxV0,  (1.8) 

which  appears  to  be  more  general  but  in  fact  corresponds  to  an  implicit  transformation  of  the 
variables  a,  0  to  another  pair  of  such  variables. 

In  terms  of  the  description  (1.6),  one  finds  that  the  nonlinear  vector  equation  (1.1)  reduces  to 
the  pair  of  scalar  nonlinear  equations 


j  •  Va  =  0,  j-V0  =  0,  (1.9) 

which  are  clearly  consistent  with  equations  (1.7)  for  a  force-free  field,  for  which  j  is  parallel  to  B. 
In  terms  of  a  and  0,  equations  (1.9)  become 

\V  x  (Va  X  V0)\  •  Va  =  0, 

(1.10) 

(V  x  (Va  X  V/))]  •  V0  =  0, 

which  look  rather  formidable  when  written  out  explicitly.  Nevertheless,  these  equations  may  be 
solved  by  a  relaxation  procedure,  which  was  the  approach  adopted  by  Sturrock  and  Woodbury 
(1967)  for  a  problem  in  rectangular  Cartesian  coordinates  and  by  Barnes  and  Sturrock  (1972)  for  a 
problem  in  cylindrical  coordinates.  In  the  simplest  version  of  the  relaxation  procedure,  one  attempts 
to  adjust  a  and  0  at  each  mesh  point  in  turn  to  arrange  that  equations  (1.10)  are  satisfied  at  that 
mesh  point.  An  alternative  procedure  is  to  attempt  to  relax  simultaneously  an  array  of  points,  such 
as  those  on  a  line,  but  this  requires  inversion  of  a  matrix.  It  was  by  no  means  clear  what  was  the 
best  numerical  procedure  for  ensuring  both  numerical  stability  and  the  most  rapid  convergence. 

Sakurai  (1979)  has  proposed  a  different  procedure  aimed  at  the  calculation  of  general  three- 
dimensional  (3D)  force-free  magnetic-field  configurations.  Sakurai  also  describes  the  field  in  terms 
of  Clebsch  variables,  but  he  uses  the  fact  that  the  force-free  field  satisfies  the  variation  principle 

0[-J-  [  d0x  (Va  x  V0)a]  =  0.  (1.11) 

8T  Jy 

Sakurai  takes  advantage  of  the  variation  principle  to  employ  the  Rayleigh-Ritz  method  (Courant 
and  Hilbert  1953)  of  representing  the  field  in  terms  of  a  set  of  base  functions  and  solving  for  the 
coefficients  arising  in  this  representation. 

In  another  article,  Sakurai  (1981)  has  discussed  another  method  for  computing  force-free  fields. 
If  fl„,  the  normal  value  of  the  magnetic  field  is  given  on  a  bounding  surface,  and  if  A  is  specified 

3 


either  in  the  region  of  the  surface  where  Bn  is  positive  or  in  the  region  where  Bn  is  negative,  the 
magnetic  field  can  be  found  by  an  iterative  process.  Starting  from  the  potential  field  with  the  given 
B„,  the  value  of  A  is  used  to  distribute  currents  along  the  field  lines.  The  field  is  then  recalculated, 
and  the  procedure  repeated  until  the  calculation  converges. 


II.  THE  MAGNETO-FRICTIONAL  METHOD 


We  have  recently  begun  to  explore  what  first  appeared  to  be  an  entirely  different  approach  to 
this  problem.  This  technique,  we  find,  is  closely  related  to  a  method  developed  by  Chodura  and 
Schlilter  (1981)  for  application  to  MHD  problems. 

Suppose  that  there  is  a  frictional  force  acting  on  the  plasma  if  it  moves  with  a  non-zero  velocity 
with  respect  to  a  fictitious  fixed  medium.  Then  the  equation  of  motion  of  the  plasma  is 


7  +«>•  Vt?)  =  -Vp  +  pg  +  P-vff, 
at 


(2.1) 


where  P  again  represents  the  Lorentz  force, p  and  p  are  density  and  pressure, &  is  the  gravitational 
field,  and  v  is  the  coefficient  of  friction.  In  the  situations  that  lead  to  force-free  fields,  p  and  p  are 
negligible,  so  that  equation  (2.1)  leads  to 


t?  =  i ~lP. 


(2.2) 


We  are  considering  the  case  that  the  plasma  has  high  electrical  conductivity  so  that 

dB  - 

—  =  Vx{vxB). 

This  equation  is  satisfied  if  we  require  that  a  and  f)  "move  with  the  plasma,”  i.e.,  that 

~  +  if-  Va  =  0, 

%  +  <r-VH  =  o. 

On  combining  equations  (2.2)  and  (2.4),  we  obtain  the  equations 

£-• 


(2.3) 


(2.4) 


(2.5) 


where  u  is  so  far  an  arbitrary  function  of  space  and  time.  In  practice,  it  has  proved  convenient  to 
choose  the  expression  v  so  that  equations  (2.5)  lead  to  the  following  increments  for  a  and  $  for  each 
time  step: 

P  ■  Va 


6a  =  — p 
6fi=-(t 


£3  ’ 
P- V0 

B *  ’ 


(2.6) 


mvmmnnamra 


where  /i  is  an  adjustable  parameter. 

The  techniques  represented  by  equations  (1.10)  and  (2.6)  for  computing  nonlinear  force-free 
fields  appear  to  be  quite  different,  yet  this  proves  not  to  be  the  case.  If  we  consider,  for  simplicity,  a 
rectangular  coordinate  system  with  mesh  size  h,  Woodbury  (1973)  showed  that  the  variations  6a, 6 f) 
require  to  satisfy  equations  (1.10)  are  given  by 

|V/9|*  -Va-Vn  (l 
Va  ■  V 0  |Va|2  J  \l 

We  rind  that  the  determinant  of  the  matrix  is  given 

A  =  |Vor|2|V/J|*  -  (1 

so  that  equation  (2.7)  leads  to 

$«\  =  h3  f  |Ve»|* 

Mj  6B*\-Vo-V^ 

This  equation  may  be  rewritten  as 

tfa\  _  h 2 

6fl)  ~  6 £» 

which  is  clearly  of  the  same  form  as  equation  (2.6). 

In  practice,  it  is  found  advantageous  to  ”under-relax”  by  introducing  into  the  coefficient  of 
(2.10)  an  additional  multiplying  constant  less  than  unity.  Hence,  although  the  relaxation  approach 
and  the  magneto-frictional  approach  begin  from  quite  different  premises,  they  lead  eventually  to 
similar  equations  for  computing  force-free  fields. 

III.  ENERGY  THEOREM 

According  to  our  model,  the  only  form  of  stored  energy  is  magnetic,  and  energy  is  continually 
decreasing  since  the  magnetic  field  lines  are  continually  doing  work  in  moving  against  friction.  If  W 
is  the  total  magnetic  energy, 

<31> 

we  see  that 

<3-2> 

On  using  equation  (2.2),  this  becomes 

<3-3> 

This  confirms  our  conjecture  that  the  magnetic  energy  will  decrease  monotonically,  and  that 
dWIdt  =  0  only  when  f  —  0  everywhere,  that  is,  the  magnetic  field  is  force-free. 


Suppose  that  the  force-free  magnetic  field,  for  given  boundary  conditions  restricting  a  and 
is  Bq-  Consider  a  perturbation  from  the  force-free  state,  the  magnitude  of  the  perturbation  being 
determined  by  the  small  parameter  w.  Then  we  may  write 


B  =  Bq  +  ZjBi  +  +  . . . 


(3.4) 


As  a  result  of  the  variational  theorem,  equation  (1.11), 


0,... 


so  that,  in  the  expansion 


W  =!*'o  +  wtt,i+w1K'2  +  ...  , 


W i=0  and  W3 

However,  since  P0  =  0,  we  see  from  (3.3)  that 


dW , 
dt 


J 


<?x 


Hence,  to  lowest  order  in  the  perturbation, 


(3.5) 


(3.6) 

(3.7) 


(3.8) 


AW  dt^  A.JdixB2' 

where 

W=W0  +  AW.  (3.10) 

It  is  clear  from  equation  (3.9)  that  those  perturbations  of  the  force-free  field  that  have  the  highest 
wave  number  will  be  those  that  damp  most  rapidly.  Hence  ’’local”  fluctuations  will  damp  rapidly, 
but  ’’global”  perturbations  will  damp  slowly. 


In  order  to  avoid  numerical  instability,  one  must  adopt  smaller  time  steps  as  one  chooses  a 
smaller  spatial  grid  size.  For  these  reasons,  it  is  advantageous  to  begin  with  a  coarse  mesh  since 
this  makes  it  possible  to  get  rid  of  perturbations  of  large  wave  number  in  a  reasonable  number  of 
time  steps.  One  may  then  choose  successively  smaller  grid  sizes  since  one  is  then  getting  rid  of 
perturbations  of  increasingly  higher  wave  number.  One  of  the  advantages  of  the  magneto-frictional 
viewpoint  is  that  it  suggests  that  one  adopt  a  regridding  procedure,  and  this  is  found  to  be  highly 
advantageous. 

One  may  specify  an  eigenvalue  problem  by  looking  for  perturbations  that,  to  lowest  order, 
satisfy  the  equations 


di  +  ffct\  —  0, 


I'rY 


1 


and  requiring  that  the  perturbations  a i,/Jj  are  zero  on  the  boundary  S,on  which  a  and  f)  are 
specified.  We  see  from  equation  (2.5)  that  these  equations  may  be  written  explicitly  as 

--^[(V  x  (Va0  x  V^o))  X  (Va0  x  +  Vax  x  V/90) 

+  (V  x  (Va0  x  V^i  +  Vori  x  V^o))  x  (Vqq  x  V^o)[  ■  Vo©  +  =  0,  (3.12a) 


-—((V  x  (Va0  x  V^0))  x  (Vor0  x  V0X  +  Vax  x  V 0O) 

+  (V  x  (Vo0  x  Vjdx  +  Varj  x  V/9o))  x  (Va<>  x  V$o)]  •  V/?o  +  <r& 1  =  0.  (3.126) 

It  is  clear  that  the  eigenfunction  corresponding  to  the  smallest  value  of  the  eigenvalue  a  is  the 
one  that  will  dominate  asymptotically  in  determining  the  rate  of  decrease  of  the  magnetic  energy. 
Hence  we  expect  that,  asymptotically,  the  magnetic  energy  will  decrease  exponentially. 

IV.  TEST  CASE 

In  order  to  study  the  variation,  during  the  iteration  process,  of  the  departure  of  the  calculated 
magnetic  field  from  the  force-free  field,  it  was  convenient  to  base  a  test  case  on  a  known  force-free 
configuration.  The  following  functional  form  for  a  and  0 

a(ar,  z)  =  ~~-e~lzcoskx, 

"  /  M  «  \ 


0(*|  V.  *)  =  V  -  7(»,  *)  =  V  -  {jj  -  1)1/2* 


yields  the  configuration 


Bz  =  --rBoe  lzco»kx, 
k 

By  =  -(1  -  ^j)1/tB0e~uco»kx, 

Bz  =  BoC~lz  tinkx, 

which  is  readily  shown  to  be  force-free.  We  adopted  the  correct  values  of  a  and  7  on  the  boundaries 
of  a  rectangle,  but  we  added  a  large  perturbation  in  the  interior  so  that  the  initial  functions  were  in 
error  by  about  a  factor  of  two. 

In  Figure  1  we  show  (a)  the  initial  trial  form  for  the  potential  a,  and  (b,c,d)  the  form  after  5,  15, 
and  60  iterations,  respectively.  Figure  1(d)  is  very  close  to  the  exact  solution.  After  100  iterations, 
the  calculated  contours  are  indistinguishable  from  the  exact  contours,  and  the  Lorentz  force  differs 
from  zero  by  only  about  10~*  of  its  initial  magnitude.  These  calculations  have  also  been  carried  out 
using  a  32  x  32  grid.  Even  in  this  case,  the  calculations  take  only  a  few  seconds  on  a  VAX/750. 

Figure  2  is  a  plot  of  the  variation  of  the  excess  of  the  estimated  magnetic  energy  over  its  exact 
value  (given  by  the  analytical  solution)  as  the  iteration  proceeds.  These  early  iterations  were  carried 
out  with  a  coarse  (8  x  8)  mesh.  We  see  that  the  convergence  is  indeed  close  to  exponential. 
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V.  CYLINDRICAL  MODEL 


Barnes  and  Sturrock  (1972)  published  the  results  of  a  calculation  of  a  cylindrically  symmetric 
nonlinear  force-free  field  configuration.  The  model  was  chosen  to  resemble  the  field  that  would  be 
developed  if  a  sunspot  were  to  rotate  with  respect  to  the  surrounding  region  of  opposite  magnetic 
polarity.  Since  the  model  adopted  by  Barnes  and  Sturrock  has  a  mild  singularity  at  the  origin,  we 
modify  it  slightly  to  avoid  that  singularity.  Using  z,  r,<j>  coordinates,  we  assume  that  a  =  0(2,  r) 
and  that 

^  =  0-7(2,  r).  (5.1) 

The  function  7  traces  the  variation  of  <fi  along  a  fi°ld  line  since  =  const,  on  a  field  line.  We  find 


Ida  „  _dadi  da 

B‘~rd 7'  Br~  rTz'  B*-Trd-z~TzTr'  (5’2) 

so  that  2 ira(z,r)  is  the  magnetic  flux  embraced  by  the  circle  2  =  const., r  =  const.,  provided  that 
we  adopt  0(2,  0)  =  0. 

We  have  adopted  the  following  boundary  conditions  for  a  and  7  on  the  plane  2  =  0: 
o(0,  r)  =  B0r2e~r’/R\ 


,  ,  /  0.5' 

3’r)  =  t-o. 


5~iMCOs(irr/2R),  0  <  r  <  2R; 


r  >  2 R. 


As  we  see  from  Figure  3,  the  magnetic  polarity  reverses  at  r  =  R,  and  there  is  a  relative  rotation 
between  the  negative  polarity  region  and  the  positive  polarity  region,  of  maximum  value  7*/. 


We  have  calculated  a  sequence  of  magnetic-field  configurations  for  the  values  1m  =  0,  0.5t,  t, 
1.5jt,  2 r,  2.5x,  and  3x.  It  was  necessary  to  close  the  region  of  calculation  by  setting  up  supplementary 
boundaries  at  z  =  L  and  r  =  L.  The  results  of  our  calculations  are  shown  in  Figure  4.  The 
appropriate  choice  of  the  outer  boundary  depends  sensitively  on  the  value  of  7 m-  The  choice 
L  =  5R  is  quite  adequate  for  7m  =  0,  but  quite  unacceptable  for  7 m  =  0.5t.  The  choice 
L  ~  20/?,  adopted  in  calculations  leading  to  Figure  4,  is  adequate  for  7m  =  2jt,  but  severely 
distorts  the  outer  field  lines  for  the  case  7m  =  3x. 


We  have  computed  the  total  magnetic  energy  for  each  case.  The  result  is  listed  in  Table  I, 
together  with  our  estimate  of  the  total  magnetic  energy  of  the  corresponding  open-field  configuration. 
We  see  that,  as  the  stress  of  the  force-free  field  increases,  the  total  energy  approaches  but  does  not 
exceed  that  of  the  corresponding  open-field  configuration.  There  is  an  important  difference  between 
these  results  of  the  present  article  and  those  of  an  earlier  article  by  Barnes  and  Sturrock  (1972). 
According  to  the  earlier  results,  it  appeared  that  the  magnetic  energy  of  a  force-free  field  could 
exceed  the  energy  of  the  corresponding  open  configuration.  Aly  (1984)  has  conjectured  that  the 
energy  of  any  force-free  field  cannot  exceed  that  of  the  corresponding  open-field  configuration,  and 
suggests  that  the  results  of  the  Barnes-Sturrock  calculations  were  in  error.  The  present  calculations 
are  consistent  with  Aly’s  conjecture.  When  one  can  be  confident  that  the  calculation  of  the  field 


configuration  and  the  field  energy  are  not  significantly  affected  by  the  external  boundary,  we  find  that 
the  energy  of  the  force-free  field  is  indeed  less  than  that  of  the  corresponding  open-field  configuration. 
Aly  in  fact  suggested  that  the  earlier  results  could  be  attributed  to  the  influence  of  the  external 
boundaries  in  the  Barnes-Sturrock  calculation. 

It  appears,  therefore,  that  our  present  calculations  do  not  support  the  proposal  of  Barnes  and 
Sturrock  (1972)  that  the  energy  of  the  force-free  field  can  exceed  that  of  the  corresponding  open-field 
configuration.  On  the  contrary,  our  present  studies  support  the  conjecture  of  Aly  (1984)  that  the 
energy  of  a  force-free  field  cannot  exceed  that  of  the  corresponding  open  field.  We  intend  to  discuss 
this  issue  in  more  detail  in  a  later  article. 

VI.  DISCUSSION 

The  magneto-frictional  method  for  computing  nonlinear  force-free  magnetic-field  configurations 
is  simple  in  concept,  fairly  simple  in  operation,  and  appears  to  offer  a  promising  approach  for  the 
calculation  of  general  force-free  magnetic-field  configurations. 

There  are  two  types  of  problems  for  which  one  needs  the  capability  of  computing  force-free 
fields.  One  of  these  is  the  calculation  of  the  field  patterns  that  will  be  produced  by  certain  models 
of  astrophysical  systems.  Our  discussion  of  the  magnetic  field  of  a  rotated  sunspot  is  an  example  of 
such  a  problem.  For  this  category  of  problems,  the  specification  of  the  field  in  terms  of  the  boundary 
conditions  of  the  Clebsch  variables  is  appropriate  and  convenient. 

However,  another  class  of  problem  is  that  of  attempting  to  model  the  force-free  field  of  a  solar 
active  region,  in  which  case  one  attempts  to  find  a  field  pattern  that  matches  certain  observational 
data.  A  vector  magnetograph  is  now  in  operation  at  Marshall  Space  Flight  Center  (  Krall  et.al. 
1982  )  that  yields  the  vector  field  of  an  active  region  at  the  photosphere.  In  order  to  compute  the 
corresponding  force-free  field,  using  the  magneto-frictional  method,  it  will  be  necessary  to  develop 
an  iteration  procedure.  First  assume  a  set  of  boundary  conditions  for  a  and  0  that  generate  the 
correct  value  of  the  known  magnetic  field  and  compute  the  corresponding  force-free  field.  Then 
compare  the  observed  transverse  component  of  magnetic  field  at  the  photosphere  with  that  which 
arises  in  the  calculation.  Based  on  that  comparison,  modify  the  boundary  conditions  for  a  and  0. 
Since  one  is  already  using  an  iteration  procedure,  and  proceeding  in  stages  by  adopting  progressively 
finer  grid,  it  is  clear  that  one  should  amalgamate  both  iteration  procedures  so  that  gross  errors  in 
the  boundary  conditions  are  corrected  at  the  coarse-grid  ievel,  and  small-scale  errors  are  corrected 
at  the  fine-grid  level. 

The  calculation  of  real  magnetic-field  configurations  also  requires  that  the  magneto-frictional 
method  be  executed  in  three  dimensions  rather  than  only  two  dimensions  as  in  this  article.  This 
will  no  doubt  require  a  more  powerful  computer  than  the  VAX/750  that  was  used  for  the  present 
calculations.  Apart  from  this  requirement,  application  of  the  procedure  to  three  dimensions  appears 
straightforward.  It  does  seem  likely  that  it  will  be  more  convenient  to  describe  the  field  in  the 


form  of  equation  (1.8)  rather  than  in  the  form  of  equation  (1.6).  For  instance,  one  could  adopt 
any  convenient  specification  of  the  variables  a  and  f),  in  either  the  positive  polarity  region  or  the 
negative  polarity  region,  and  then  adjust  the  quantity  f  to  match  the  observed  value  of  the  normal 
magnetic  field. 

In  addition  to  the  extension  of  this  procedure  to  three  dimensions,  the  general  approach  seems 
appropriate  for  application  to  more  general  field  configurations  than  that  of  a  simple  force-free  field. 
For  instance,  in  solar  active  regions,  it  often  happens  that  one  or  more  current  sheets  develop.  It 
should  be  possible  to  divide  the  region  occupied  by  the  field  into  two  regions,  in  each  of  which  the 
magnetic  field  is  force  free,  but  allowing  for  a  possible  current  sheet  at  the  interface  between  the  two 
regions.  In  principle,  one  should  be  able  to  apply  the  magneto-frictional  method  to  each  of  these 
regions  separately.  However,  the  interface  must  be  treated  as  a  free  boundary.  If,  at  each  point  on 
the  boundary,  the  magnetic  field  strength  is  different  on  the  two  sides  of  the  boundary,  then  the 
boundary  should  be  moved  in  response  to  the  net  pressure  acting  on  the  boundary  at  that  point. 
By  moving  the  boundary  progressively  in  accordance  with  this  rule,  one  should  tend  asymptotically 
to  the  state  in  which  each  region  is  force  free  and  the  interface  has  the  appropriate  equilibrium 
configuration.  We  hope  to  explore  some  of  these  and  other  developments  in  future  articles. 

A  shortcoming  of  the  present  computer  program  is  that  one  must  calculate  the  field  over  a  very 
large  region  when  the  field  is  highly  stressed.  This  is  due  to  the  tendency  of  the  force-free  field  to 
inflate,  moving  towards  an  open-field  configuration.  A  possible  solution  to  this  problem  is  to  adopt 
a  larger  mesh  size  in  the  weak-field  region.  Another  possibility,  which  would  be  well  suited  for  a 
wide  range  of  problems,  is  to  adopt  a  spherical  coordinate  system  and  arrange  that  the  increment 
in  radius  be  proportional  to  the  radius.  We  shall  report  on  the  results  of  this  approach  in  a  later 
article. 

This  work  was  supported  in  part  by  the  Office  of  Naval  Research  Contract  NOOOH-85-K-Olll, 
NASA  Grants  NGL05-020-272  and  NAGW-92. 
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TABLE  1 

MAGNETIC  FIELD  ENERGY  AS  A 
FUNCTION  OF  ROTATION 


1.5x 

2r 

2.5x 

openfield 


(1)  This  estimate  is  only  approximate  due  to  slow  convergence  and  proximity  of  the  boundary. 
The  estimate  for  =  3x  is  even  less  certain  and  is  not  given. 


(2)  The  energy  of  the  open  field  has  been  estimated  by  adding  to  the  energy  of  the  field  within 
the  box  an  estimate  of  the  energy  of  the  field  outside  the  box. 
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Figure  Captions 

Figure  1.  Contours  of  constant  a  (field  lines)  for  the  two-dimensional  problem  described  in 
Section  IV.  (a)  is  the  assumed  initial  form,  and  (b),  (c),  and  (d)  give  the  forms  after  5,  15,  and  60 
iterations,  reaspectively.  These  calculations  were  carried  out  on  a  very  coarse  (8  x  8)  grid. 

Figure  2.  The  excess  energy  of  the  computed  force-free  field  over  the  energy  of  the  exact  field, 
as  a  function  of  iteration  number,  for  the  case  desribed  in  Section  IV. 

Figure  3.  Contours  of  constant  a  for  the  problem  of  cylindrical  symmetry  described  in  Section 
V.  (a),  (b),  (c),  (d),  (e),  (f),  and  (g)  give  the  contours  for  the  following  values  of  7m,  the  relative 
rotation  of  the  regions  of  positive  and  negative  polarity:  0,  0.5jt,  x,  1.5x,  2j r,  2.5x,  and  3x.  Figure  3(h) 
gives  the  contours  of  constant  a  for  the  corresonding  open-field  configuration. 

Figure  4.  The  energy  of  the  cylindrical  force-free  field  model  of  Section  V  as  a  function  of  7m, 
which  measures  the  relative  rotation  of  the  positive  polarity  region  and  the  negative  polarity  region. 
Also  shown  is  the  energy  of  the  corresponding  open-field  configuration. 
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